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Abstract We present a classification of non-liermitian random matrices based on 
implementing commuting discrete symmetries. It contains 43 classes. 
This generalizes the classification of hermitian random matrices due to 
Altland-Zirnbauer and it also extends the Ginibre ensembles of non- 
hermitian matrices [1]. 

Random matrix theory originates from the work of Wigner and Dyson 
on random hamiltonians [2]. Since then it has been apphed to a large 
variety of problems ranging from enumerative topology, combinatorics, 
to localization phenomena, fluctuating surfaces, integrable or chaotic 
systems, etc... Non- hermitian random matrices also have applications 
to interesting quantum problems such as open choatic scattering, dissi- 
pative quantum maps, non-hermitian localization, etc... See e.g. ref. [6] 
for an introduction. The aim of this short note is to extend the Dyson 
[2] and Altland-Zirnbauer [4] classifications of hermitian random matrix 
ensembles to the non-hermitian ones. 

1. What are the rules? 

As usual, random matrix ensembles are constructed by selecting classes 
of matrices with specified properties under discrete symmetries [2, 3]. To 
define these ensembles we have to specify (i) what are the discrete sym- 
metries, (ii) what are the equivalence relations among the matrices, and 
(iii) what are the probablility measures for each class. 
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(i) What are these discrete symmetries. 
Let h denote a complex matrix. We demand that the transformations 

specifying random matrix classes arc involutions — their actions are of 
order two. So we consider the following set of symmetries: 



C sym. : 


h = 




(1) 


P sym. : 


h = 


~php^^, = 1 


(2) 


Q sym. : 


h = 


q q^^, q^q~^ = 1 


(3) 


K sym. : 


h = 


kh*k-^, kk* = ±l 


(4) 



denotes the transposed matrix of h, h* its complex conjugate and 
its hermitian conjugate. The factor Cc is just a sign ec = ±. We could 
have introduced similar signs in the definitions of type Q and type K 
symmetries; however they can be removed by redefining h — > ih. 

We also demand that these transformations are implemented by uni- 
tary transformations: 

cc^ = 1, pp^ = 1, qq^ = 1, kk^ = 1 (5) 

In the case of hermitian matrices one refers to type C symmetries 
as particle/hole symmetries or time reversal symmetries depending on 
whether ec = — or Cc = + respectively. Matrices with type P symmetry 
are said to be chiral. Both type Q and type K symmetries impose reality 
conditions on h and they are redundant for hermitian matrices. 

(ii) What are the equivalence relations. 

We consider matrices up to unitary changes of basis, 

h ^ uhu^ (6) 

In other words, matrices linked by unitary similarity transformations 
are said to be gauge equivalent. For the symmetries (1-4), this gauge 
equivalence translates into: 

c —> ucv^ , p^upu^^, q ~^ uqu\ k —> uku*~^ (7) 

The classification relies heavily on this rule and on the assumed uni- 
tary implementations of the discrete symmetries. 

We shall only classify minimal classes, which by definition are those 
whose matrices do not commute with a fixed matrix. 

(iii) What are the probability measures. 

Since each of the classes we shall describe below is a subset of the space 
of complex matrices, the simplest probability measure ii{dh) one may 
choose is obtained by restriction of the gaussian one defined by 

^l{dh) = M exp(-Tr^/i^ ) dh (8) 
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with M a normalization factor. It is invariant under the map (6). 

There is of course some degree of arbitrariness in formulating these 
rules, in particular concerning the choice of the gauge equivalence (6). 
It however originates on one hand by requiring the gaussian measure 
(8) be invariant, and one the another hand from considering auxiliary 
hermitian matrices Ti, obtained by doubling the vector spaces on which 
the matrices h are acting. These doubled matrices are defined by: 

o) (^) 

They are always chiral as they anticommute with 75 = diag(l, —1). Any 
similarity transformations h uhu~^ are mapped into 7i — > UTCU^ with 
U = diag(tt, ti^ ^^). So, demanding that these transformations also act 
by similarity on Ti, imposes u to be unitary. 

On Ti, both type P and Q symmetries act as chiral transforma- 
tions, n ^ -rnV'^ with V = diag(P,P) and H QHQ'^ with 

Q = ^ ^ ^ ) ' ^ may be block diagonalized if is Q or P sym- 
metric. Indeed, if /i is Q symmetric then Q and H. may be simultaneously 
diagonalized since they commute. If /i is P symmetric, H. commutes with 
the product 75^. 

Type C and K symmetries both act as particle/hole symmetries re- 
lating Ti to its transpose Ti.^ . The classification of the doubled hamilto- 
nians H thus reduces to that of chiral random matrices, cf. [4]. However, 
the spectra of h and Ti may differ significantly so that we need a finer 
classification involving h per se. 

2. Intrinsic definition of classes. 

To specify classes we demand that the matrices belonging to a given 
class be invariant under one or more of the symmetries (1-4). It is im- 
portant to bear in mind that when imposing two or more symmetries it 
is the group generated by these symmetries which is meaningful. Indeed, 
these groups may be presented in various ways depending on which gen- 
erators one picks. For instance, if a matrix possesses both a type P and 
a type C symmetry, then it automatically has another type C symmetry 
with d = pc and = —Ec- 

The intrinsic classification concerns the classification of the symmetry 
groups generated by the transformations (1-4). 

We demand, as usual, that the transformations (1-4) commute. For 
any pair of symmetries the commutativity conditions read: 

c = zizpcp^ ; p* = zkk~^pk ; q = zizpqp^ (10) 
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The signs it are arbitrary; they shall correspond to different groups. 
Without reality conditions. 

This arises if no type P and type Q symmetry is imposed so that 
no reality condition is specified and h is simply a complex matrix. We 
may then impose cither a type P or a type C symmetry or both. Not 
all groups generated by a type P and a type C symmetry arc distinct 
since, as mentioned above, the product of these symmetries is another 
type C symmetry but with an opposite sign Sc- The list of inequivalent 
symmetry groups, together with the inequivalent choices of the sign ec, 
is the following: 



Generators 


Discrete symmetry group 


Number 




Defining relations 


of classes 


No sym 


no condition 


1 


P sym 


p^ = l 


1 


C sym 


= ±c, ec 


4 


P, C sym 


= 1, = ±c, pcp^ = c 


2 


P, C sym 


p^ = 1, = ±ec c, pcp^ = — c 


2 



If the sign e^, does not appear as an entry it means that the value of this 
sign is irrelevent — opposite values correspond to identical groups. The 
sign factors =b written explicitly are relevant — meaning that e.g. the 
groups generated by a type C symmetry with = c ox = —c are 
inequivalent. The cqtiivalcnces among the defining relations for groups 
generated a type P and a type C symmetry arc the following: 

{p^ = 1,C^ = ±C,C = pcp'^)e^ = {p^ = 1,C^ = ±C,C = pcp'^)-ec 

(p^ = 1,0^ = ±ec c, c = -pcp'^)ec — {p^ = l,c^ = ±ec c, c = -pcp'^)-e^ 

In the above table we anticipate the numbers of classes which we shall 
describe in the following section. They depend on the numbers of in- 
equivalent representations of each set of defining relations. 

Considering discrete groups generated by more symmetries of type P 
or C does not lead to new minimal classes. For instance, if the group is 
generated by two type P symmetries, then their product commutes with 
h and thus they do not define a minimal class. Similarly, suppose that we 
impose two type C symmetries with sign ed and €c2- If ^ci^c2 = —, their 
product makes a type P symmetry and the group they generate is among 
the ones listed above. If eci^c2 = +> their product commutes with h and 
they do not specify a minimal class. More generally, considering more 
combinations of type P and type C symmetries does not lead to new 
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minimal classes as in such cases one may always define fixed matrices 
commuting with h. 

With reality conditions. 

This arises if we impose at least a type Q or a type K symmetry, but 
we may simultaneously also impose symmetries of other types. Again, 
there could be different but equivalent presentations of the same discrete 
group as not all of these symmetries are independent. For instance, a 
product of a type P symmetry with a type C, Q oi K symmetry is again 
a type C, Q or K. The list of inequivalent groups generated by two or 
three of these symmetries is the following: 



Generators 


Discrete symmetry groups 


Number 




Defining relations 


of classes 


Q sym 


q = q^ 


2 


K sym 


kk* = ±1 


2 


P, Q sym 


= 1, = 1, q = ±pqp^ 


2 


P, K sym 


/ = 1, kk* = ±1, kp* =pk 


2 


P, K sym 


p^ — 1, kk* = 1, kp* = —pk 


1 


Q, C sym 


q = q\ = ±c, q^ = c'q^'^c, tc 


8 


Q, C sym 


q = q\ = ±c, = -c+9"^c, ec 


4 


P, Q, C sym 


= 1, q = q\ = ±c, ec 






c = tcppcp^ , q = epqpqp\ q^ = €cq c^q~^c 


12 



As above, the explicitly mentioned signs it correspond to inequivalent 
groups. The groups with defining relations (p^ = 1, kk* = 1, kp* = 
—pk) or (jp = 1, kk* = —1, kp* = —pk) are equivalent. The groups 
generated either by a type Q and a type K, or by a type C and a type 
K symmetries are included in this list because the symmetries of type 
C, Q or K are linked by the fact the product of two of them produces 
a symmetry of the third type. The last cases, quoted in the last line of 
the above list, are made of groups generated by three symmetries, one 
of a type P and two of type either C, Q or K. Their defining relations 
depend on the choices of the signs Ccp, ^pq and Ccg. The equivalences 
between these choices are the following: 

(c — €cp — £pq — £cq — ~l~)ec ~ — ^cp — ^PQ — ^cq — i)— 

(c — ^cp — ^P<1 — ^c<l — ~l~)ec ~ (/' — -^^1 ^cp — ^P<1 — ^c<l — i)— 

(c — ic; tcp — ~^pq — ^cq — i)ec — ~ ~^cp — ^pq — ^cq — ~)— Cc 

(c — i€c C; 6cp — ^pq — ^cq — i)€c — ('^ ~ i^c €cp — Cpq — Ccg — )~€c 

(c — iCc C; €cp — ^pq — — )ec — ~ i^c C; Ccp — tpq = €cq — )— €c 

Considering groups generated by more symmetries does not lead to new 
minimal classes since in such cases one may construct matrices commut- 
ing with h. 
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3. 



Explicit realizations of the classes. 



Having determined the inequivalent groups of commuting discrete 
symmetries, the second step consists in finding all inequivalent repre- 
sentations of the defining relations of those groups. Due to the rules we 
choose, in particular the second one, eq.(6), we only consider represen- 
tations in which all symmetries are unitarily implemented and which are 
not unitarily equivalent. 

Wc shall list all these representations, adding an index ec to recall 
when the action (1) of the corresponding discrete group depends on Ec- 

Without reality conditions: 

If we impose no symmetry at all, the class is simply the set of complex 
matrices. 

If we impose only a type P symmetry. The matrix p is unitary and 
square to the identity, so that it is diagonalizable with eigenvalues ±1. 
The solution p = 1 is trivial as it implies ^ = 0. Assuming for simplicity 
that the numbers of +1 and —1 eigenvalues are equal, we may choose a 
basis diagonalizing p: 



Here and below, o"^, (j^ and ay denote the standard Pauli matrices. 

If we only impose a type C symmetry, c is unitary and either sym- 
metric and antisymmetric. As is well known [2], up to an appropriate 
gauge choice it may be presented into one of the following forms: 



Let us now impose a type P and a type C symmetry. In the basis 
diagonalizing p with p = tr^ 1, the commutativity relation pcp^ = ibc 
means that either p and c commute or anticommute. If they commute, 
then c is block diagonal in this basis so that c=l(8)lorc=l(8' icTy 
depending whether it is symmetric or antisymmetric. If they anticom- 
mute, c is block off-diagonal so that, modulo unitary change of basis, it 
may be reduced to c = cxa; (8" 1 or c = iuy ® 1. However, as explained in 
the previous intrinsic classification, these two cases correspond the same 
symmetry group but with opposite signs Cc- So a set of inequivalent 
representations is: 





(c=l)6c; (c = ^o■y(8)l)e,; 



(12) 



(p = (g) 1, c = 1 ® 1); {p 



fj^®!®!, c=l(8) iUy 1); 




(13) 



Altogether there are 10 classes of non-hermitian random matrices with- 
out reality conditions. They are of course parallel to the 10 classes of 
hermitian random matrices [4]. 
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With reality conditions: 

Imposing only a type Q symmetry, we have q = and qq^ = 1 since, 
by choice, we implement the symmetry unitarily. So q is diagonalizable 
with eigenvalues ±1. The solution g = 1 is non-trivial as it simply means 
that h is hermitian. When q possesses both +1 and —1 eigenvalues we 
assume that they are in equal numbers, hence: 

(g = 1); (g = a, 1); (14) 

Imposing only a type K symmetry, we have kk* = ±1, kk'^ = 1, and 
their classification is similar to that of type C symmetries. There are 
two cases: 

(k = 1); {k = iay «) 1); (15) 

When imposing both type P and Q symmetries, p and g both square 
to the identity and either commute and anticommute. In the gauge in 
which p = az ® 1 the possible g are 10 1, cr^ ^ 1 or cr^ (g) 1, up to 
unitary similarity transformations preserving the form of p. However, 
the two first possibilities generate identical groups, thus the inequivalent 
representations are: 

(p = (g) 1, g = 1 (g) 1); (p = (8) 1, g = cTj; (g) 1); (16) 

Let us next impose a type P and a type K symmetry. With the unitar- 
ity property of p and k, the commutativity relations between these sym- 
metries are similar to those between type P and C symmetries. Thus, 
the classification of these representations may be borrowed from eq.(13) 
and we have: 

(p = (g) 1, A; = 1 (g 1); (p = (72 <g 1 <g 1, A; = 1 (g iay 1); 

(p = (72(g)l, A; = (7^; (g) 1); (17) 

When imposing simultaneously a type Q and a type C symmetry, 
their product is a symmetry of type K. The commutativity constraints 
between type Q and C symmetries are solved in the same way as the 
commutativity conditions for type P and C symmetries; the only differ- 
ence being that g = 1 is now a non-trivial solution, contrary to p = 1. 
The inequivalent representations for g and c are: 

(g = l, c=l),, ; (g = lg)l, c = i(7y®l)e,; (18) 

(g = (72 (g 1, C = 1 (g l)e^ ; (g = (72 (g 1, C = 1 (g i(Ty)e^ 

(g = (72 <g 1, c = iay (g l)g^ ; (g = (72 (g 1, c = cr^, (g l)g^; 
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Finally, let us impose a type P symmetry together with two among the 
three types of symmetries Q, C and K. The commutativity constraints 
are solved by first choosing a gauge in which p = az 1- The list of 
inequivalent solutions is: 

{p = az I, g = l(8)l, c = l(X'lorc=l(8) io'y); 
{p = az®l, q=l(g)l, c = «) (19) 
= (T2 (8> 1, Q = (Ta; (8> 1, c = 1 (8) 1 or c = 1 (g) iay 

01 c = ax <^ 1 or c = ax ic^y)ec] 

Here, each choice of c corresponds to a different class. 

Altogether, eqs. (11-19) give 43 classes of random non-hermitian ma- 
trices. As for the Ginibre ensembles [5] — which correspond to the 
classes without symmetry or with type K symmetry with A; = 1 or 
k = iay (8> 1 — we expect that in the large matrix size limit the density 
of states in each class covers a bounded domain of the complex plane 
with the topology of a disc. 
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